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§1.  Introduction 

Being  among  the  simplest  models  in  the  kinetic  theory  of  gases,  the  Broadwell  model  has  served 
as  a  paradigm  to  understand  the  phenomenon  of  relaxation  and  the  transition  from  a  microscopic 
to  a  macroscopic  description  of  gases.  It  consists  of  the  system  of  semilinear  hyperbolic  equations 

and  derives  from  a  six-velocity  model  when  specializing  to  one-dimensional  flows,  for  which  the 
densities  of  particles  moving  in  directions  orthogonal  to  the  flow  are  all  equal.  (We  refer  to 
Broadwell  [B]  or  to  Platkowski  and  Diner  [PI]  for  the  derivation  in  the  kinetic  theory  context). 
The  function  /  =  (/i,/2,/3)  is  defined  for  (x,i)  €  E.  X  E'*’  and  describes  densities  of  particles:  /i 
for  particles  moving  in  the  positive  ^-direction,  /z  in  the  negative  ^-direction  and  /j  in  each  of  the 
positive  or  negative  p-  or  ar-directions;  as  a  consequence  we  confine  to  solutions  with  fj  >  0.  The 
parameter  e  stands  for  the  mean  free  path,  a  measure  of  the  average  distance  between  successive 
collisions.  The  right  hand  side  in  (1.1)  is  called  the  collision  operator  and  measures  the  rate  of  gain 
(or  loss)  in  densities  of  particles  effected  through  collisions.  It  is  characterized  by  the  quantity 

Q(/)  =  /|-/i/2.  (1.2) 

The  zeroes  of  Q(f)  are  the  states  of  equilibrium  for  the  system,  /I  =  /1/2,  and  are  called 
Maxwellians.  Finally,  associated  with  each  /  are  the  quantities  p/  =  /i  +  f 2  +  4/3,  m/  =  /i  —  /2, 
measuring  the  local  density  and  momentum  flux  in  the  a;— direction,  respectively. 

The  limit  when  the  mean  free  path  approaches  zero  is  known  as  the  fluid  dynamic  limit.  For 
small  mean  free  path  the  strong  interactions  of  particles  allow  a  macroscopic  description  of  the 
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flow  to  become  meaningful.  In  the  case  of  the  BroadweU  model  the  induced  macroscopic  ’’Euler 
equations”  are  easy  to  identify.  First  rewrite  (1.1)  in  the  form 

^(/i  +  /j  +  4/3)  +  ^(/i  -  h)-^i 

^(/l-/2)+^(/l+/2)  =  0,  (1.3) 

9/3  ^  f  2  f  f  \ 

It  is  formally  expected  that  cts  e  0  the  first  two  equations  remain  unchanged  while  the  third 
causes  the  limit  of  /  to  be  a  local  Maxwellian.  If  we  denote  by  Ihe  limit  of  /,  then 

and  F  =  (jFi,F2)  satisfies  the  limiting  fluid  equations 


+  F2+  4{FiF2y^^)  +  ^{Fi  -  F2)  =  0 , 


The  corresponding  macroscopic  density  and  momentum  of  the  fluid  are  given  by 

p  =  +  4(FiF2)^/2  +  F2  ,  m  =  pu  =  -  F2  . 


The  algebraic  system  can  be  easily  inverted  and  leads  to  an  alternative  form  of  the  limit  ’’Euler 
equations”,  in  terms  of  the  macroscopic  variables  (p  ,  u). 


If + 


where  g{u)  :=  5  p(l  +  -  ij .  They  form  a  strictly  hyperbolic,  genuinely  nonlinear  system 

of  conservation  laws  (Caflisch  [C]). 

The  justification  of  the  fluid  dynamic  limit  has  been  the  object  of  several  investigations.  We 
refer  to  Cercignani  [Ce]  for  a  survey  of  the  literature  on  the  Boltzmann  equation  and  to  Platkowski 
and  Diner  [PI]  for  results  on  discrete  velocity  models  of  kinetic  theory.  For  the  BroadweU  model, 
the  fluid  dynamic  limit  is  understood  for  smooth  solutions  of  the  limit  ”Euler  equations”  (Inoue 
and  Nishida  [NI],  Caflisch  and  Papanicolaou  [CP]).  Regarding  the  case  of  solutions  with  shocks,  we 
mention  the  studies  on  stability  (in  time)  of  traveUng  wave  solutions  (Kawashima  and  Matsumura 
[KM])  or  rarefaction  wave  solutions  (Matsumura  [M])  for  the  BroadweU  model,  and  a  recent  study 
by  Xin  [X],  showing  that  a  given  piecewise  smooth  solution  with  noninteracting  shocks  of  the  Umit 
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fluid  equations  can  be  approximated  by  solutions  of  the  Broadwell  system  as  e  ^  0,  that  gives 
a  definitive  answer  to  one  direction  of  the  problem.  The  converse  problem,  to  show  that  a  given 
family  of  solutions  to  the  Broadwell  system  converges  globally  in  time  to  a  fluid-dynamical  solution, 
remains  at  present  open. 

Insight  in  the  latter  direction  is  provided  by  the  approach  of  self-similar  fluid  dynamic  limits, 
introduced  in  Slemrod  and  Tzavaras  [ST]  and  further  studied  in  Tzavaras  [T]  and  Fan  [Fj.  For 
Riemann,  Maxwellian  data 

/(i,0)  =  {  {_’  *  ^  ^  with  (1-7) 

if,  X  <0 

W-)  =  W+)  =  o,  >0,  (M) 


admissible  solutions  of  the  limit  fluid  equations  (1.4)  are  expected  to  be  self-similar  functions  of 
^  =  x/f.  On  the  other  hand,  the  Broadwell  system  does  not  possess  space- time  dilational  invariance 
and  does  not  admit  self-similar  solutions  of  that  type.  Motivated  by  an  analogous  idea  for  systems 
of  conservation  laws  (Dafermos  [D]),  one  considers  a  modified  Broadwell  system 


9fi  ,9fi  _  1  2  f  f^. 

^  -  2A  -  L(  f'i  -  f  f  \ 
dt  dx 


(1.8) 


which  does  preserve  the  invariance  under  dilations  (x,  t)  —>■  (ax,  at),  a  >  0.  The  problem  (1.7  — 1.8) 
admits  self-similar  solutions  of  the  form  /  =  f{x/t),  that  are  constructed  by  solving  the  singular 
boundary  value  problem 


=  for  {€[-1,1) 

-a,  =  -^(fl  -  fih) 

A-i)  =  r ,  /(+!)  =  /* 

subject  to  data  satisfying  (M).  There  are  three  goals  to  be  attained:  (i)  To  construct  solutions 
of  the  problem  (Pg)  for  e  >  0  fixed,  (ii)  To  show  that  as  £  ^  0-f  the  solutions  converge 
a.e  to  a  local  Maxwellian  F  solution  of  the  Riemann  problem  (1.4,  1.7).  (iii)  To  investigate  the 
structure  induced  on  F  from  emerging  through  self-similar  fluid  dynamic  limits. 
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In  this  article  we  present  a  survey  of  applications  of  the  method  to  the  Broadwell  model  ([ST], 
[T]  and  [F])  and  comparisons  with  studies  of  shock  profiles  ([B],  [C]).  The  hope  is  the  approach 
will  be  useful  in  studying  relevant  questions  for  other  kinetic  theory  models  or  for  hyperbolic 
conservation  laws  with  relaxation  terms.  The  presence  of  relaxation  mechanisms  is  natural  in 
many  physical  contexts  and  has  been  investigated  extensively  {e,g.  Liu  [L],  Chen,  Levermore  and 
Liu  [CLL]),  both  in  the  realm  of  kinetic  theories  as  well  as  for  models  that  arise  in  various  branches 
of  continuum  physics.  One  objective  is  to  obtain  a  quantitative  understanding  of  the  regularizing 
effect  induced  on  shocks  by  relaxation.  In  this  direction,  we  mention  the  comparison  of  shock 
profiles  with  traveling  wave  solutions  of  an  associated  system  of  viscous  conservation  laws  arising 
from  (1.1)  via  the  Chapman- Enskog  expansion  [Cj.  An  important  difference  of  the  self-similar 
relaxation  investigated  here  is  that  it  penalizes  the  whole  wave  fan  simultaneously.  Comparisons 
between  the  structure  entailed  by  self-similar  limits  and  the  Broadwell  shock  profiles  are  carried 
out  in  the  text. 

§2.  The  limiting  fluid  equations 

First  certain  properties  of  the  limit  fluid  equations  (1.4)  for  solutions  Fi  >  0  are  presented.  These 
properties  are  discussed  in  Caflisch  [Cj;  we  give  an  independent  presentation  for  completeness  and 
to  account  for  extended  differences  in  notation.  The  characteristic  speeds  Ai(F),  ^2(F)  are  the 
solutions  of  the  binomial 


(Fi  +Fj  +  -  (fi  -  f!)Ai  -  v/JTfi  =  0; 


(2.1) 


they  are  real  and  are  given  by  the  expressions 

1 


F,-F2±  \l(F,  +  Fi+2^/Fjif  -4Fji 


2  (Fi  +  ^2  +  \/FiF2) 
where  Ai  corresponds  to  the  minus  sign  and  A2  to  the  plus.  One  easily  checks 

Fi  -  F2 


-1<Ai(F)<0<A2(F)<  1,  Ai(F)< 


Fi  -j-  F2  4\/ Fi  F2 


The  corresponding  right  eigenvectors  are  given  by 


n  = 


A,-  +  1 
Ai-1 


i  =  1,  2. 


(2.2) 


=  u  <  A2(F)  (2.3) 
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Upon  differentiating  (2.1)  we  can  express  in  terms  of  A*,  and  use  the  resulting  relations 

to  compute 


VXi-ri  = 


(1  -  A.-2)  (  (Fi  +F2+  4v/j\Jl)A.  -(Fi-F2) 


)u.  ( 

.F2  \ 


)■ 


(2.4) 


2^/KF^  \iFi  +  F2+VF[r\)2Xi-{Fi-F2), 

It  follows  from  (2.2),  (2.3)  and  (2.4)  that  VAi  •  r,-  >  0  for  i  =  1,  2,  and  the  system  (1.4)  is  strictly 
hyperbolic  and  genuinely  nonlinear. 

The  shock  curves  for  (1.4)  are  defined  by  solving  the  Rankine-Hugoniot  conditions,  expressed 
cis  the  equivalent  algebraic  system 


-s  ( «i  +  2  «3  )  +  ui  =  -s  ( tni  +  2 103  )  +  ini 

-s  (  tl2  +  2  W3  )  -  U2  =  -«  (  102  +  2  IO3  )  -  102  (2.5) 

—  U1U2  =  -  lOi  IO2  =  0 

for  the  states  u  =  (iti,ii25M3)  and  to  =  (101,102,103)  and  the  shock  speed  s.  To  this  end  fix  one 
state,  say  lo,  with  103  =  y'101102  and  consider  the  increments  [itj]  =  Uf  —  lo^  Then 

2  s  [113]  =  (!-«)  [tti]  =  -(1  +  s)  [112]  (2.6) 

[«3]  ([tts]  +  2 103)  =  [«i]  [112]  +  102  [lii]  +  wi  [U2]  (2.7) 


Substitution  of  (2.6)  into  (2.7)  yields 


[«3] 


(1  +  3s^)  [113]  -  2  [(loi  +  102  +  103)5^  -  (loi  -  102)5  -  103 


=  0 


If  [113]  =  0  there  is  no  shock.  Thus,  using  (2.1)  and  (2.6),  we  obtain  a  representation  of  the  shock 
curve  paxametrized  by  the  shock  speed  5  in  the  form 


l«3l  =  (.  -  A.(»))  (.  -  A,(»)) . 


25 


25 


(2.8) 


=  [«2]  =  -^[n3]. 

Given  a  triplet  (u,  10;  5)  satisfying  (2.8)  there  are  two  associated  shock  solutions  with  speed  5  of 
(1.4)  one  with  left  state  10  and  right  state  u  and  one  with  the  states  reversed. 


§3.  Shock  profiles  for  the  Broadwell  system 

The  regularizing  effect  induced  on  shocks  by  relaxation  can  be  understood  by  studying  traveling 

wave  solutions  for  (1.1)  (c.f.  Broadwell  [B],  Caflisch  [C]).  If  we  introduce  the  ansatz 

,x  —  st.  X  -  st  /  .  X  ,+  ^  X 

fj  =  ^  ’  i;j(±oo)  =  ff  (3.1) 
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and  look  for  a  traveling  wave  connecting  two  end  states  /-  and  then  vj  satisfy  the  system  of 
differential  equations 


-(s+i)^  =  Q(v), 


(3.2) 


After  some  simple  manipulations,  it  turns  out  that  v  =  (wi,  V2>  ^3)  must  satisfy  the  set  of  equations 


2s  {v3{t)  -  /3  )  =  (1  -  s)  (ui(r)  -  )  =  -(1  +  s)  (v2{t)  -  /a  ) , 


dV3 

dr 


2s 


(ul  -  V1V2), 


(3.3) 

(3.4) 


that  u(±(xi)  must  be  Maxwellian  states  satisfying  the  Rankine-Hugoniot  conditions  (2.5),  and  thus 
the  jumps  [fj]  —  —  fj  of  the  data  are  connected  with  the  shock  speed  s  through  relations  (2.8). 

Equation  (3.4)  may  be  written  in  the  form 


dvz 

dr 


l  +  3s^ 
2s(l  -  s^) 


(1^3  -  h  )(^3  -  fz), 


(3.5) 


the  solution  of  which  is  given  by  the  explicit  formula 

«3(r)  =  ta-nh  -  fa  )  (’’  “  ’'o))  (3.6) 

where  the  constant  tq  determines  the  shift  of  the  shock  profile.  Then  ui(r)  and  V2{t)  are  obtained 
from  (3.3). 

The  question  arises  which  shock  solutions  of  (1.4)  have  associated  shock  profiles.  From  (3.6),  a 
triplet  (/■*",/“;  s)  satisfying  (2.8)  admits  a  shock  profile  with  f~  as  left  state  and  /■*■  as  right  state 
if  and  only  if  s(.s^  -  l)(/3’  -  f^)  >  0.  It  follows,  upon  using  (2.8)  to  express  [fz]  (with  w  = 
and  with  w  =  f^)  and  (2.3),  that  the  latter  condition  is  equivalent  to 


Ai(/‘'’)  <  s  <  Xi{f  )  for  s  <  0 , 
A2(/‘'’)  <  s  <  A2(/“)  for  s  >  0 . 


(3.7) 


These  are  the  Lax  shock  conditions  for  (1.4). 


§4.  Self-similar  fluid  dynamic  limits 

We  turn  now  to  our  main  objective  the  study  of  the  system  (Ve)  and  its  limits  £  — >■  0.  Regarding 
the  question  of  existence  of  solutions  we  have 


7 


Theorem  1  (Slemrod  and  Tzavar<is  [ST],  Fan  [F]).  For  any  satisfying  (Af)  and  e  >  0  the 
boundary-value  problem  {Vg)  has  a  solution  continuously  differentiable  in  (— 1,0)U(0, 1)  and 
Holder  continuous  with  exponent  0  <  <  1  of  the  singular  points  ^  =  — 1,0,+1.  The  regularity 

improves  as  e  decreases  and  for  e  <  So  the  function  is  Lipshitz  continuous  on  [-1, 1]. 

The  functions  are  extended  to  the  whole  real  line  by  setting  /®  =  f~  on  (-oo, -1)  and 
ff  =  on  (l,oo).  Regarding  the  self-similar  fluid-dynamic  limit  we  have 

Theorem  2  (Slemrod  and  Tzavaras  [ST]).  Let  {/*}e>o  he  a  family  of  extended  solutions  of  (Ve) 
corresponding  to  data  subject  to  (M).  Then: 

(i)  There  are  positive  constants  mj ,  Mj  and  Kj  ,  j  =  1, 2, 3,  depending  on  the  boundary  data 
but  independent  of  e  such  that 

^€[-1,1]  (4.1) 

//  <  Ki  (4.2) 

(ii)  There  exists  a  subsequence  {ff"}  with  0  and  a  positive,  bounded  function  F  of  bounded 

variation  such  that  f^"  —y  F  pointwise  on  the  reals.  The  function  F  satisfies 

Fs  =  y/ Fi  F2  for  a.e.  ^  €  [—1,1]  (4.3) 

and  the  balance  of  mass  and  momentum  equations 

-^■^AFt  -I-  ^2  -h 4(^-1  F2)^/2)  -h  ^(Fi  -  F2)  =  0, 

-f-(F’i-F’2)-H-(Fi-l-F2)  =  0, 

m  the  sense  of  measures. 

Remarkably  the  problem  of  existence  is  more  difficult  than  performing  the  £  — >  0  limit.  There 
are  two  approaches  to  tackle  the  existence  question.  In  [ST]  a  Fredholm  alternative  type  of  the¬ 
ory  for  singular  boundary  value  problems  is  developed.  The  theory  has  a  wide  range  of  potential 
applicability,  but  at  the  final  stage  one  must  show  that  a  certain  linear  but  non- autonomous 
boundary-value  problem  has  no  eigensolutions.  As  there  are  no  general  techniques  for  such  ques¬ 
tions  this  is  a  technical  obstacle.  For  the  Carleman  model  this  method  works  for  any  Maxwellian 
data,  but  for  the  BroadweU  restrictions  on  /*,  of  technical  nature,  had  to  be  imposed.  These  were 
removed  in  [F],  where  the  existence  question  is  addressed  by  first  desingularizing  the  system  and 
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then  taking  a  dynamical  systems  point  of  view  for  the  resulting  boundary-value  problem.  Fan  [F] 
fixes  and  uses  a  shooting  method  to  construct  solutions  on  (—1,0)  and  (0,1)  separately,  and 
then  shows  that  the  traces  of  the  constructed  solutions  at  ^  =  0  intersect  transversally. 

The  uniform  variation  estimates  are  based  on  the  following  observation.  The  collision  operator 
Q{f^)  satisfies  the  differential  equation 


dQin n 

^  ^  +  1 


(4.5) 


The  uniqueness  theorem  for  differential  equations  implies  that  on  each  of  the  intervals  (—1,0)  and 
(0, 1)  either  Q(/®)  is  identically  zero  or  it  never  vanishes.  {Ve)  in  turn  implies  that  the  /J  are 
monotone  or  constant  on  the  respective  intervals,  and  (4.1),  (4.2)  foUow  by  a  case  analysis  and  a 
use  of  the  balance  of  mass  equation  for  the  bounds.  One  also  obtains  the  uniform  estimate 

f\9in\ir<C  (4.6) 


HeUy’s  theorem  implies  there  exists  a  subsequence  and  a  function  F  of  bounded  variation 

so  that  “*>  F  pointwise  on  IR.  Use  of  (4.1)  and  (4.6)  shows  that  the  components  Fj  >  0  and  F 
is  a  local  Maxwellian.  Passing  to  the  limit  e  0  in  {Vc)^  we  deduce  F  satisfies  (4.4)  in  the  sense 
of  measures. 

For  a  function  of  bounded  variation  the  right  and  left  limits  F(^~),  F{^+)  exist  at  any  ^  and 
its  domain  can  be  decomposed  into  two  disjoint  subsets  :  C  the  points  of  continuity  of  F,  and  S 
the  points  of  jump  discontinuity  of  F,  with  S  at  most  countable.  The  components  Fj  inherit  the 
monotonicity  properties  of  /J.  On  account  of  (4.6) 


F3(o  =  {m)m)f\  eec, 

F3(^±)=  (f’i(^±)F2(^±))'/\ 

(4.4)  implies  that  at  any  point  ^  G  <S  the  Rankine-Hugoniot  conditions  are  satisfied 


-^([Fi  -l-Fa-f  |^!)+  [Fi  -F2]  =  0 

-H[Fi-F2]\ll)+[Fr  +  F,]\ll  =  0 


(4.8) 


By  construction  F  =  /_  on  (-00, -1]  and  F  =  /+  on  [l,oo).  The  function  ^Fi(|),F2(f)^  is 
easily  seen  to  be  a  weak  solution  of  the  Riemann  problem  (1.4),  (1.7). 
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§5.  Structure  induced  by  fluid  dynamic  limits 


In  the  previous  section  we  outlined  the  construction  of  a  solution  to  the  Riemann  problem  (1-4), 
(1.7)  via  self-similar  fluid  dynamic  limits.  In  this  section  we  discuss  the  structure  of  the  resulting 
hmit,  emphasizing  the  behavior  at  points  of  discontinuity.  We  outline  the  main  ingredients  of  the 
technique  and  refer  to  Tzavaras  [T]  for  details  of  the  presentation. 

First,  the  appropriate  framework  for  passing  to  the  f  — ^  0  limit  in  is  that  of  measures. 
Consider  the  functions 

$^(0= 

J  —  OO  ^ 

and  note  that  takes  constant  values  outside  [-1, 1].  The  sequence  is  of  uniformly  bounded 
variation  on  IR.  Helly’s  selection  principle  implies  there  exists  a  subsequence,  denoted  again  by 
and  a  function  of  bounded  variation  $  such  that  $  pointwise  on  IR.  In  turn,  Helly’s 

convergence  theorem  implies 


(5.1) 


<  v,  ip  >  : 


J 


J 


(pd^  =:<  v,  p> 


(5.2) 


for  any  tp  £  Cc(lR-),  continuous  function  with  compact  support.  By  the  Riesz  representation 
theorem  u‘,  v  may  be  viewed  as  finite  (signed)  Borel-Stieltjes  measures,  both  supported  on  [—1, 1], 
is  generated  by  and  v  is  generated  by  $,  the  right  continuous  version  of  $  defined  by 
$(x)  =  $(x+).  Equation  (5.2)  states  that  v‘  v  va.  the  weak-*  topology  of  measures.  It  allows 
to  pass  to  the  limit  £  -+  0  in  {Ve)  and  obtain,  for  any  test  function  (p  £  C'g(lEl), 


J  Fi((^-l)yj)'d^=<i/,V5> 

j  F2  ((^  +  l)¥>)'df  =<v,<p>  (5.3) 

jF^{i^)'d^  =-i  <u,^> 

with  Fz  =  y/FiFz.  It  is  suggested  by  (5.3)  and  can  be  justified  by  an  analysis  near  the  singular 
points  that  supp  u  is  precisely  the  set  of  points  where  /  is  not  a  constant  state.  Equation  (5.3) 
implies  the  weak  form  (4.4),  but  carries  additional  information  if  further  properties  of  u  are  known. 

The  second  ingredient  is  a  representation  formula  for  Q{f^)le.  Using  (4.5)  and  (Ve)  it  is  shown 
in  [T]  that 

QiD  _  I  (/i  (0)  -  /r )  Mi  on  (- 1 , 0) 

(fiW- on  (0,1) 


e 


(5.4) 
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where  //i,  fj,^  are  given  by 


/^i(0  = 


/°1  ei^‘(0  dC  /",  ^  exp{i  c^(s)  ds}  dC 

1  rC 


-C 

eis‘(0 


li  c+T  Si  C+T  ®^P^  e  /i+ 


^€(-1,0), 


ee(o,i), 


(5.5) 

(5.6) 


and  a_,  a+  are  any  fixed  points  with  — 1  <  a_  <  0  <  a+  <  1.  Note  that  because  of  their  form 
(along  subsequences)  ^  /x^.  and  //I  — *  fi-  weak-*  in  measures.  We  deduce  from  (5.2)  and  (5.4) 
that  i>  arises  as  a  limit  of  probability  measures  and  supp  v  C  supp  /ii_  U  supp  fj,^. 

Remarkably,  the  function  c®  in  (5.5  —  5.6)  and  (4.5)  is  connected  to  the  wave  speeds  of  the 
hyperbolic  system  (1.4),  as  it  can  be  expressed  in  the  form 


(1-a^ 


(/f + /I + /I)  -  Mn)  (e  -  Hn) , 


(5.7) 


where  Xi^2{f)  ^^re  given  by 
Ai,2(/)  = 


fl-f2±  \/(/i+/2+ 2/3)2  -4/1/2 


(5.8) 


2  (/i  +  /2  +  /a)  L 

A  comparison  of  (5.8)  and  (2.2)  shows  that  Xi,2if)  coincide  with  the  wave  speeds  Ai,2(i^)  along 
Maxwellian  states.  Hence,  along  the  convergent  sequence  /®  — F, 

1 


c  = 


(Fi  +F2  +  Fa)  (^  -  Ai(F))  (e  -  A2(F)) 


(5.9) 


(1-f^)^ 

pointwise  on  (—1,0)  U  (0, 1). 

Use  of  the  above  ingredients,  together  with  an  analysis  of  the  behavior  near  the  singular  points 
gives  ([T],  Lemma  4.3,  Proposition  4.4) 


Proposition  3  .  There  are  constants  Ai±,  A2±  with  —1  <  Ai_  <  Ai+  <  0  <  A2-  <  A2+  <  1  such 
that  suppi/  C  (supp/x_  U  supp^.j.)  C  [Ai_,Ai+]  U  [A2_,A2+].  Moreover 

(i)  If  ^  ^  supp  /x+  then  c(s)ds  <  c(s)ds  for  any  (  €  (0, 1) 

(ii)  If  i  &  supp/i_  then  c{s)ds  <  c{s)ds  for  any  (  €  (—1,0). 


The  maximization  properties  (i)  and  (ii)  capture  the  effect  induced  on  shocks  by  self-similar 
fluid  dynamic  limits.  To  illustrate,  fix  ^  <  0,  ^  €  «S.  Then  ^  €  supp/x_,  the  function  g  achieves  its 
global  maximum  in  (-1,0)  at  ^  and  thus  c(^-|-)  <  0,  c(^— )  >  0.  Using  (5.9)  and  (2.3)  we  deduce 
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^i{Fi^+))  <  ^  a  weak  form  of  the  Lax  shock  conditions.  In  fact  the  complete 

behavior  of  F  can  de  characterized  and  the  final  result  is  stated  in  the  first  part  of  Theorem  4. 

The  last  topic  is  to  discuss  the  relation  between  self-similar  limits  and  shock  profiles.  Let  ^  £  S 
a  point  of  discontinuity  for  F,  and  note  that  F{^—)  5^  F{^+)  are  Maxwellian  states  satisfying  (4.8). 
Given  a  sequence  of  points  with  estimates  (4.1  —  4.2)  imply  the  functions  {v|}  defined 

by  Vj(0  =  //(^e  +  sC)  are  of  uniformly  bounded  variation  in  the  new  variable  —  00  <  ^  <  00. 
This  accounts  for  a  shift  of  the  shock  in  the  original  solution  and  the  introduction  of  the  stretched 
variable  C-  Helly’s  theorem  and  a  diagonal  argument  in  turn  imply  the  existence  of  a  subsequence 
and  a  function  v  =  {vi,V2,v^)  such  that 

fji^e  +  £  C)  ^j(C)  pointwise  for  —  00  <  (  <  00  .  (5.10) 

Then  vj  are  solutions  of  the  traveling  wave  equations  (3.2)  (c./.  second  part  of  Theorem  4). 

Theorem  4  (Tzavaras  [T]).  The  limiting  function  F  constructed  in  Theorem  3  has  the  behavior: 
F  stays  constant  on  each  connected  component  0/  IR  —  supp  u.  Also  there  exist  disjoint  closed  sets 
I\^,  with  each  h,  associated  with  the  k-th  characteristic  speedy  such  that  supp  i/  =  7;^^  U  and 

(i)  either  Ixf^  is  empty,  or 

(ii)  Ixk  contains  a  single  point  in  S,  in  which  case  F  is  a  shock  wave  on  Ix,^  satisfying  (4.8)  and 
the  Lax  shock  conditions,  or 

(iii)  Ixi,  is  a  full  interval  of  points  in  C,  in  which  case  F  is  a  k-rarefaction  wave  on  Ix,^. 

Given  a  point  of  discontinuity  f  G  there  exists  a  choice  of  the  sequence  such  that  t;(Q  in 
(5.10)  is  a  (smooth)  solution  of  (3.2)  that  satisfies  lim  t;(C)  =  F{(-),  lim  t;(^)  =  F((+). 

C— *--*00 

We  remark  that  the  proof  is  of  purely  analytical  nature.  Geometric  properties  such  as  the 
genuine  nonlinearity  of  (1.4)  and  the  form  of  the  shock  curves  (2.8)  are  only  used  in  the  last  step 
to  exclude  contact  discontinuities  and  simplify  the  emerging  structure  of  F. 
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